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ABSTRACT 

It has been argued that the solar magnetic cycle consists of two main periodic com- 
ponents: a low-frequency component (Hale's 22-year cycle) and a high-frequency com- 
ponent (quasi-biennial cycle). The existence of the double magnetic cycle on the Sun is 
confirmed using Stanford, Mount Wilson and Kitt Peak magnetograph data from 1976 
to 1996 (solar cycles 21 and 22). In the frame of the Parker's dynamo theory a model 
of the double magnetic cycle is presented. This model is based on the idea of two dy- 
namo sources separated in space. The first source of the dynamo action is located near 
the bottom of the convection zone, and the second operates near the top. The model 
is formulated in terms of two coupled systems of non-linear differential equations. It 
is demonstrated that in the case of weak interaction between the two dynamo sources 
the basic features of the double magnetic cycle such as existence of two component and 
observed temporal variations of high-frequency component can be reproduced. 

Subject headings: Sun: magnetic field — Sun: activity — Suns interior 
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1. Introduction 

The 22-yr cycle of solar activity is a magnetic cycle 
which consists of two 11-years sunspot cycles. The 
22-yr cycles begin on even sunspot cycles accord- 



ing t o the Zurich numbering ( Gncvyshcv and Ohl 



1948 ) and manifest themselves in reversal of polar- 
ity of sunspots (Hale's law). From one 11-yr cycle 
to another, the polarity of the preceding ( p ) and 
following ( / ) sunspots reverses. This reversal cor- 
responds to the changing of the toroidal magnetic 
field ( B v -component). During even- numbered cy- 
cles B v > (the direction of the magnetic vec- 
tor coincides with the direction of solar rotation) in 
the northern ( N ) hemisphere and B v < in the 
southern ( S ) hemisphere. This relation is reversed 
during odd-numbered cycles. In parallel with this, 
the poloidal magnetic field, or background field ( B r - 
component), shows a 22-yr periodicity (Howard and 
LaBonte, 198l|), and other studies flStenflo, 199d|) 



have suggested that the low-latitude and polar fields 
similarly show a 22-yr periodicity . Furthermore, 
these studies have hinted at a shorter term periodicity 
of about 2 yr (high-frequency component). In addi- 
tion, during the periods of the three-fold polar mag- 
netic field reversals which occur during some sunspot 
cycle maxima, the temporal separation of the zones of 
alternated polarity of the magnetic field o n H a charts 



is ap p roximately equal to 1. 5-2.5 years ( Waldmeier 



1973 , Benevolenskaya, 1991 ). Similar periodi cities 



were found in vari ations of radio flux on 10.7c m (Bel- 



mont et al., 1966), flares and sunspot areas (Akioka 



t al., 1987). The high-frequency, 2-year, component 



is substantially weaker than the main 22-year cycle 
and its intensity varies with time. The two year com- 
ponent more clearly appeared in northern hemisphere 
in cycle 20 and in southern hemisphere in cycle 21; its 



value was smaller in cycle 22 (Benevolenskaya, 1996). 
However, the biennial cycle represents a challenge to 
the solar dynamo theories which usually explain only 
the main cycle. In this Letter, I present a new evi- 
dence of the two components of the solar cycle, and 
an attempt to explain the results in term of Parker's 



dynamo theory (Parker, 1979) 



2. Analysis of magnetograph data 

Using Kitt Peak magnetograph data we show how 
synoptic (longitude independent) magnetic field pat- 
terns evolve through cycles 21 and 22. Each of the 
synoptic maps is represented by the values of Bu as a 



function sine latitude and Carrington longitude. The 
Carrington coordinate system is a reference frame ro- 
tating rigidly at the rate which corresponds to the 
synodic rotation rate of sunspots at a latitude of 
about 16° . The observed line-of-sight components 
( -B|| ) averaging over all longitudes for each Carring- 
ton Rotation is represented in Figure lb. The time 
step in the time set of values Bu corresponds to Car- 
rington rotation period ( P — 27.2753 days). 

The relative sunspot number is shown in Figure la. 
Component B r (Figure lc), was calculated from ob- 
servational data assuming that the true average field 
direction is radial and Bg = . The radial compo- 
nent will be easily found at all latitudes besides near 
the poles: 



B r = B||/sin0, 



(1) 



where 9 is colatitude. 

To separate the high-frequency component in the 
data, we apply a difference filter by computing AB r / At 
for different time intervals At . This is a reasonable 
procedure if the poloidal magnetic field consists of two 
components and is represented as two dynamo waves: 



B r (t,9) = b r {9)smuj c t 



A 



sin(uJ8t + ip). (2) 



where b r (9) is the amplitude of the low- frequency 
radial component of the magnetic field; A is ra- 
tio between amplitudes of low-frequency and high- 
frequency components; lj c is the frequency of the 
Hale's solar cycle; los is the frequency of biennial 
cycle; tp is a phase. 

The expression for the AB r / At can be written as 



&B r (t,9) , b r (0) 

— b r (9)Lo c cosoj c t H — lu$ cos{Lust + <p). 



At 



A 



(3) 

When los — I0ui c , A = 2 as it took place in cycle 
20, then the low- frequency term dominates in the ex- 
pression for B r (equation ||) and the high-frequency 
term prevails in the expression for AB r j At (equa- 
tion %. Thus, presuming the existence of the double 
magnetic cycle, one finds that the low-frequency term 
prevails in the B r -component, while the high - fre- 
quency term dominates in AB r j At . 

Contour plots of AB r /At as a function (or 
cos 9 ) and time, t, measured in Carrington rotation 



2 




76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 











^ — 














































pi* c \ y-w-i 




























SipWr ■ ::r 

















Fig. 1. — (a) The values of the sunspot number 
for cycles 21 and 22 as function of time (in years), 
(b) The line-of-sight component of the solar mag- 
netic field as function of Carrington rotation num- 
ber, (c) The radial component of the solar magnetic 
field (blue lines correspond to isolines for +3, +10 
G , red lines, for -3, -10 G). (d) The values of the 
AB r /At for At = 6 solar Carrington rotations (blue 
lines correspond to isolines for +0.3, +3 G /rotation, 
red lines, for -0.3, -3 G /rotation). 



are represented in Figure Id for At = 6P . Dur- 
ing the cycles 21 and 22 the zones of increasing and 
decreasing strength of the surface magnetic field ap- 
pear in both the N and S hemispheres (see Figure 
Id). The width of these zones is approximately 2 
years. This result coincides with investigations done 
earlier based on Wilcox and Mount Wilson synop- 
tic maps. Corresponding power spectra for AB r / At 
show that the component with period T = 1.5 — 2.5 
years is dom i nant in both hemisph eres ( Bencvolcn-| 
skaya, 1995 , Benevolenskaya, 1996 ). These results 



are the basis of our suggestion of a two component 
of solar cycle. 

3. Model 



The possible explanation of the double magnetic 
cycle is that the magnetic fields are generated by 
Parker's dynamo acting in convective zone. The 
low-frequency component is generated at the base of 
the convective zone due to large scale radial shear 
dQ jdr ; £1 is the angular velocity. The high-frequency 
component may be generated in subsurface regions 
due to latitudinal shear d£l/d8 or due to radial shear. 
The recent investigations of solar interior rotation 
show a significant radial gradient of angular velocity 
exists in subsurface of the convective zone together 
with the latitudina l gradient of the angular velocity 
( Bchou, ct al., 199§1 ) . For simplicity we use only latitu- 
dinal shear for generation of the high-frequency com- 
ponent. Cartesian coordinates are employed, with x 
denoting the radial, y the azimuthal and z the lat- 
itudinal coordinate. We consider axisymmetrical so- 
d 

lutions ( — = ) . 



dtp 

At the base of the convection zone turbulence is 



suppressed by strong magnetic field (Parker, 1993) 
and, therefore, diffusivity in the first layer 771 will 
could be less then r\2 in the second layer. The ax- 
isymmetrical mean magnetic field is decomposed into 
toroidal and poloidal parts. The following set of equa- 
tions describe evolution of the magnetic field in thin 
layers at the levels x\ and X2 '■ 
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dA yl 

at 

dB yl 
dt 

dCi 
~df 



dA 



y2 



dt 

dB V 2 
dt 



(ai + Ci)B yl + i]iAA yl 



^G x + m AB yl 



-vC\ +pA yl B yl 

(a 2 + C 2 )B y 2+mAA y 2 
dA y2i 



(4) 



dx 



-G z + m AB 



!)2 



(5) 



dC_ 2 
dt 



—^r = -vC 2 + pA y2 B y 2 



y is the toroidal magnetic field and A y is 



where B 

the azimuthal vector potential which gives poloidal 
field, r\ is diffusivity and a is a kinetic helicity, C 
is the variable part of kinetic helicity. The first and 
the second equations in the both systems (||) and (jj|) 
are a generation of Parker's equations for large-scale 
dV 

radial shear G x = — — and large-scale latitudinal 

dx 

dV 

shear G z — correspondingly. We have used the 
dz 

dC 

third equation in the both systems — — = — vC+pAB 

dt 

according to Kleeorin and Ruzmaikin (1982) in a sim- 
ple form, for feedback of the magnetic field on the 
helicity. In these equations v > and p < are 
parameters of the feedback. These two non-linear 
differential equation systems describe the evolution 
of two independent sources of the magnetic field. In 
the frame of this model it is difficult to explain ob- 
served variations of the high-frequency component. 
Therefore, we assume that the erupted low-frequency 
magnetic field can influence the physical conditions 
in the region where high-frequency component oper- 
ates, through modifying the helicity there. This al- 
lows the high-frequency component to vary with time. 

dC 2 



In this case, the equation for C2 becomes 



dt 



—vC<x + p(A y2 + aA y i){B y 2 + aB y i) , where the pa- 
rameter a captures the feedback of the low-frequency 
component of the magnetic field on the helicity in re- 
gions where the high-frequency component is found. 
Following Weiss et al. (1984) the solutions of these 
equations were found as A\ = a(t)e lkx , B% = b(t)e' lkx 



and G\ = c a (t) + c(t)e 2lkx , Im c Q = for x\ level 
and the same expressions for X2 level where x is 
replaced by z . A, B,C are complex functions as 
A = ai + ia 2 , B = bi + ib2 ■ 

It is convenient to present the system in dimcn- 
sionless units. Since there will be an underlying 
periodicity in all calculations, we use those periods 
as relative units of time. In relative units t — ► 

t = r]2k 2 ta , a > a = a/r]2k 2 a , a = r]l/r]2, 

C — ► C = C/rfek a 0) G z — ► G z — G z /r)2ka and 
G x — > G x — G x /r]2ka , where a is some mean 
value of the kinetic helicity. Therefore, G x and G z 
are represented as dynamo numbers in our model. 
Therefore, set of the partial differential equations re- 
duced to the following equations were then solved nu- 
merically. 

= (ai + Ci)Bt - aAi + |Bf C x 
= -iA\ G x - aB 1 ' (6) 

= -uCx+pAxBx 

= (02 + C2)B 2 - A 2 + |S| C 2 
= iA 2 G z -B 2 ' (7) 

= -vC 2 +p(A 2 + aA 1 ){B 2 + aBi) 

We have investigated solutions of these two non- 
linear systems as functions of the parameter a . Re- 
sults for three principal cases, a — , 0<a<2 and 
a = 2 are represented in Figure 3, where a\ — 012 = 
1 , v — 0.5 , p — — 1 , a = 0.1 , G x = 0.3 , G z = 20 , 
which are in the reasonable range for the Sun. If 
a = we have two independent sources of the mag- 
netic field which together generate low-frequency and 
high-frequency signatures (Fig. 2 a, b) . For a weak in- 
teraction between these sources ( < a < 0.2 ) we still 
get a low and high frequencies response. (Fig. 2 c, d). 
Moreover, the high- frequency component becomes 
more regular, and consequently is decreased when a 
increases. With further increase of the feedback of the 
low-frequency magnetic field on the helicity near the 
top surface the high-frequency component disappears 
and only a low-frequency mode regime is established 
(Fig. 2 e, f). Therefore, in the case of weak interaction 
between two sources of magnetic field it is possible to 
obtain a stable double magnetic cycle. 

4. Conclusions 

Our model simulates the double magnetic cycle 
and temporal variations of the biennial cycles from 
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one 11-year cycle to another. According to the mag- 
netograph data (Benevolenskaya, 1996) the 2-year 
component more clearly appeared in the northern 
hemisphere in cycle 20 and in the southern hemi- 
sphere in cycle 21. In cycle 22 it was present in both 
hemispheres but was of lower amplitude. 

In our model, a reduced high-frequency compo- 
nent occurs when the erupted magnetic field of the 
main (Hale's) cycle imposes the helicity in the regions 
of generation of the high-frequency component. The 
high-frequency component is more pronounced when 
the effect of erupting magnetic fields of the main cycle 
on helicity in this region is small. Thus, this simple 
model provides a qualitative explanation of the double 
magnetic cycle. The next step to a quantitative model 
is in development of a 2D model in spherical geom- 
etry, since this simple model gives only a qualitative 
picture of the two-components of the solar magnetic 
cycle. 
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Fig. 2. — (a) Real part of A = A\ + A2 as function 
time for independent systems (a = 0). (b) Power 
Spectrum for case (a) in relative units, (c) Real part 
of A — A\ + A2 as function time for a — 0.17. 

(d) Power Spectrum for case (c) in relative units. 

(e) Real part of A = A\ + A2 as function time for 
a = 0.2 . (f) Power Spectrum for case (e) in relative 
units. 
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